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ABSTRACT 



Research on wind and temperature profiles over land have 
established a number of fundamental relationships that are tested 
here over an oceanic location. Most attempts at calculating the 
roughness parameter over oceanic surfaces have been based on 
surface-layer theories, and often employ an iterative procedure in 
approaching the sea surface. This paper employs a similarity 
relationship between the wind and temperature profiles, using data 
at ship "Victor" appropriate to forced convective theory, as 
defined by Priestley. 

The theory of forced convection, as applied here, makes it 
necessary to compute the Richardson number and the Monin-Obukhov 
scale length as preliminaries to the computation of roughness 
parameter. The roughness parameter and friction velocity are 
obtained from the theory using an iterative procedure. The wind 
speed at 20 meters is correlated against scale length and roughness 
parameter. The latter makes very little contribution to the variance 
of the windspeed while the scale length contributes in a manner 
reasonably expected from theory. 
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1 . Introduction 



Studies of the roughness parameter, Z Q( have been numerous 
over land and a number of fundamental relationships using various 
meteorological parameters have been determined. Far less 
research has been accomplished in this field over the oceans. 

This is readily understandable when one considers the expense and 
difficulty in erecting and maintaining an instrumented platform 
at sea . With this in mind , one could turn to what should be a 
prime source of data — the ocean station weather ships. 

Certain problems become evident very quickly when 
considering station Victor. Four different ships occupy this 
station and only two have the same anemometer height. This, at 
first, seemed to present an insurmountable problem in using the 
data but the theoretical error involved in using a median height of 
twenty meters instead of actual heights is less than two and one 
half percent. When considering that winds are recorded only to 
the nearest knot and temperature to the nearest degree, the error 
introduced in the wind profiles by considering all heights equal to 
twenty meters is insignificant. 

Monin and Obukhov C3J , in their universal similarity theory 
introduced three scaling parameters that are considered essentially 
invariant with height for layers near the ground. These three 
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parameters are scaling velocity, u*, scaling length, L, and 
scaling temperature, T*, and have been well described in a 
recent text (Lumley and Panofsky) M on atmospheric turbulence. 
L, u* and T* are defined as follows: 



I - _ Uf C P P T~ 

L - * f/T 



( 1 ) 



~r* — ~ -r-H 

1 — ^ J J* 






(2) 



U* = 



r 



(3) 



It should be noted that a selection of cases has been limited 
to cases for which -Ri^. 032 holds. According to Priestley's dri^i- 
rion Cfj our restriction limits the study to the cases of "forced 
convection". Here Ri refers to the turbulent parameter 

Ri - HMl* m 

and is to be distinguished from the flux Richardson number Rf^ 
defined as 

Rf = p- Ri ' (5) 
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Under normal circumstances, it is considerably easier to deal 
with Ri than with Rf , since the factor is itself an 

increasing function of Ri. 

2. Important properties of the forced -convective surface layer 
According to the similarity theory for the near-neutral, or 
forced convection case (as defined in section 1), the wind 
profile can be written 

U = jf (#*% + *if) 

In a like manner the temperature profile can be written, after 
Lumley and Panc^foky 

*-a=T*(]R,£+ fr) 

In both (6) and (7), Lumley'-andcPanofafey £l] suggest the value 
$ = 4.5. 

By dividing equation (7) by equation (6)^ and taking the 
derivative of both numerator and denominator with respect to 
height, one has a relationship governing the constants T* and 

uVk : 



T* _ aa/az 

W iu/ ^ 



Relationship (8) means that the change in the temperature 
over a given height is proportional to the change in the wind over 
the same height. Even more important, this ratio is constant 
with height since T , u , and k are all considered constant in the 
surface layer. 

The ratio of these constants T *A */k) may be determined 
from the air temperature at psychrometer height (six meters) and 
the temperature at Z Q , assumed to be the same as the sea 
surface temperature, so that equation (8) may be written, in 
finite-difference form centered at three meters as 




Here u Q is the wind speed at Z D and is zero. By equation (8) this 
ratio is constant at all levels in the surface layer, including that 
at 1.41 meters. 

From the definition of L, it follows that L satisfies the 
identity 




However, Friestley^4jfhas introduced a surface layer scale 
lenght L' given by 




(ID 
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so that (10) becomes 



tf = R' s 



(12) 



a formula lacking the troublesome ratio which has been 

found to depend upon Ri. But from the definition of T*, we have 



From (13), together with the similarity profiles (6) and (7), which 
apply in the case of forced convection, it follows that 
and therefore L=L' . This is a fundamental property characteristic 
of similarity profiles in forced convection. 

3. Derivation of a specific value for (Ri)j 41 

From equations (4) and (9), the Richardson number is 
expressible in the form 



where every parameter in (14) is now known except for the 




or 




(13) 




(14) 
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The surface-layer wind profile is assumed to be valid to 
heights of twenty meters. The windspeed ug may then be related 
to U 20 by a new parameter through 



where Kg_ 2 Q is formed from the ratio of Ug to U 2 Q. This now 
leaves the quantity §Z )j 41 as so ^ e unknown in 
equation (14). One may approximate this quantity by forming the 
ratio of two wind differences centered at 1.41 meters. In the 
following equations, subscripts denote the height of a variable 
in meters although substitution in centimeters is generally 
required. For example we have 




(15) 




(16) 



and 




(17) 



Therefore, it follows that 




(18) 



Now the last two unknowns in equation (14) are 'found in terms of 
known or measured quantities. Following Martin Pi 




(19) 
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and, similarly, 




with all values in centimeters. Kg_ 2 Q may be approximated 
satisfactorily for obtaining Ri by omitting In Z 0 . 

The right sides cf(l5) and (16) are now known as a function of L 
so that (14) may be solved for Ri as 




?z /a -a) 

&o U2.0 



( 21 ) 



where Kj 4 , and Kg_ 2 g are given as functions of L in (19), (20). 

In the forced convection regime, Lumley and Panofs(?y show 
that the Richardson number may be written as 

z 

Ri — I I / 8 Z ^ — (22) 

' ' H L 

Eliminating Ri between (21) and (22) leads to a quadratic in L, 
the only meaningful root of which is given by 

L - ~ - ^ - (23) 

A = \ WLm)-.7U%- c iooo Ri'iLzmLzooo) 

B = Z% Ri'(JUZOO)(l~ZOOO) (24) 

C = HUlRiY(JU2MJL ZOODY8WO ‘) 

-W^Ri'(LZOO)(L Z00O)J 



where 
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Ri may now be obtained by solving equation (22) with Z=141.1 
cm and the value of L obtained from (23) . 

The quadratic equation defined by (23) and (24) had as 
input data only windspeeds at 20 meters, and temperature - 
differences measured from the surface to 6 m. Only one data-set 
was selected at ship Victor per day, and these were selected with 
a view toward insuring cases of forced convection: for example, 
a vertical temperature difference less than or equal to and 

windspeeds as large as possible under these circumstances were 
selected. A table of Ri values computed as a function of both U 20 
and is included in the Appendix. 
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4. Iterative procedure of solution of roughness parameter 



Following Martin C. 2 J, use is made of equation (6) to obtain 
the quantity U 2 g - Uq ^ in the form 

= i f(z-3ozq2.3o\o^)+^s( 3 ^°X i £) ( 25 ) 

But with U 2 Q as follows 

U» = 2.302 + 

one obtains the result, U G j, given in (27) 

U,, = .005 U 20 + if(2. 2 l&XI-jtjftjs Zo ) ( 22 ) 



By using the Monin-Obukhov normalized wind shear 

-Avh u 

- U*-2- ,3/ 



and (12), it follows that 




Because Z Q is small compared to 20 meters, the quantity 
( c>U/<^J^r iZ^ will be rewritten involving a larger spam, as 

follows 

Ui-u , 

UmX,' Z h 



9 



where 4 will be defined immediately below. 

(£££z ) M - - Cl.H (Uzo-Uo.1) 



From (6) , we thert-have 



( j i 64 \? 

Cl.M - - LtfjMi 

Combining (28) and (29) leads to 



^ Uo.,) 



where is defined in (29). 

A third equation used in the iteration-process is 



(29) 



(30) 



U>.i “ Zo )(\ -h ) 

but this has already been built into equation (25). The latter may 
be used for more accurate determination of Z Q once an earlier 
iterate, or satisfactory guess, has been established. Equation 
(31) assumes the use 6f equation (6)~to an elevation of 10 cmTabove 
mean sea level. 

To adopt an efficient iterative method which does not lay 
too much emphasis on the last 10 cm, we must give proper 
weight to the layer 0-20 m. Here, we use equations (27) and (30) 
with u Q 1 eliminated. This leads to the result 
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and, of course, we also have 

U 20 Z ^ 200 -j- j 

The iteration process may be summarized as follows: 

1. Assume Z Q ^ = 0.1 and obtain a corresponding estimate 
of (u*/k)^ from (32). 

2. With this first estimate of (u*/k)^, use equation (33), 

(21 

with U 20 and L known, to determine Z Q V ' . 

3. Repeat the procedure using Z 0 ^ in (32) to get a second 
estimate of (u*/k)^ 2 ^ . 

4. Repeat the iteration until a final . residue is obtained: 

| z^-zH & \C>' 

When this has been accomplished Z Q is taken to be .Z and 
is the appropriate value from either (32) or (33). 

One safeguard must be followed in the use of (32): the 
denominator of the right side must not be allowed to become zero 
or negative. However, the use of (33) helps to restrict (u*/k) to 



reasonable, non-negative values. 



The results, giving Z Q , u*/k as well as U 2 q» Ri' , L and Rl 
are listed in Table 1 (Appendix). In addition, some statistical 
relationships between the parameters are also considered in the 
following section. The raw data considered was drawn from 30 
randomly selected dates. The value of L is computed from 
equation (23), while Rij ^j is obtained from (22). 

The computations were made using data cards obtained from 
the National Weather Records Center, Asheville, North Carolina, 
for ocean station "Victor". All cards selected were for 0300Z, the 
observation nearest noon local time, from the months of November 
and December of 1960 and of March 1961. The number of 
intervals between successive days was considered a random 
function of the date. Wind values were tabulated in this paper in 
centimeters per second, and temperatures in degrees? Celsius. Sea 
conditions during the period were not considered, with the 
exception that no major storms were in the immediate area. 

5 . Re s u 1 tesancb cort c 1 usAq i o n s 

The values of U 2 Q/ L and Z Q for the sample of 30 cases were 
subjected to a multiple regression analysis leading to a result of 
form 



li tt =• -30-MZo, - 0.0053 \Lt7S7.V <34> 

variable 3 variable 1 variable 2 
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Of the two independent variables tested in (34), Z 0 plays an 
insignificant part in that its partial regression coefficient was 
-0.00969. When the mean, value of Z G = 0.1066 cm was inserted 
into (34), the simple regression equation 

Uzo = - 0.00 531 L +754.76 

resulted. The variable L had a partial correlation coefficient of 
-0.56153, so that its effect upon the reduction of the variance 
of U 20 was sufficiently strong that resulting F-value due to the 
regression involving L was 6.576. This F-value, with 2 and 27 
degrees of freedom, is significant at a level in excess of the 99 
percent confidence level. This result is not surprising from the 



manner in which l/L enters U 20 E^33. A brief 
the other statistics obtained are listed: 


summary of some of 


l-U 


- 18707 cns 


Z.-OAOkbZ cn 


I =-H 3,773.6c* 




= 3 781 7c 7is 


H£ o = 0.l0 2J8cn 


<JT = H0,i82.Lcn 






15,.*= -0.00461 


1 7^-/ - ~ 0- I5~ 3 
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The most surprising fact regarding the statistics, and 
equation (34) is the complete unimportance of the variable Z Q in 
describing U2o* At first glance, equation (33) would suggest a 
much stronger dependence, since -lnZ Q is one of the terms of the 
right side, and this should indicate that U 20 decreases with 
increasing Z Q . This is a result presented by many writers, who 
emphasize the tendency for u to decrease with increasing surface 
Rossby number Vg/fZ 0 . However, from equation (33) it can be 
seen that 









LL 



1 



(36) 



where (u*/k) n is the neutral value of (u*/k). That part of (u*/k) 
which depends upon stability is primarily incorporated in the ratio 






so that the neutral part of (u*/k) is essentially 



m/c» 



The final point to be emphasized here then is as follows: 
Within equation (36) the term -lnZ Q tends to reduce U2g- However, 
carefully measured results of Hay (1955) indicate that 

(lf%= 1 3. Ig z 0 
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so that u may be expected to increase in the boundary layer as Z Q 
increases. Obviously these two effects oppose one another, with 
a wind-increasing tendency from (u*/k) and a wind-decreasing 
tendency arising from -lnZ 0 within the bracket of (36). This 
oddity resulting from opposite effects produced by an increase in 
Z Q appears to explain the lack of correlation of U 20 with Z Q . 

No statistical analysis involving u*/k with Z Q has been made 
but their values are listed in the Appendix^and may serve as the 
basis of future research on the subject. 
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APPENDIX I 



Table of selected micrometeorological parameters in the surface layer. 



If. 


AT(°C) 


Ri' (10 


8 ) L 


Ri 


u*/k 


Z o 


1081.5 


-3.88 


-1.15 


-25197.329 


-.0057 


112.510 


.098703 


412.0 


-4.44 


-9,02 


-4900.552 


-.0327 


56.065 


.250295 


309.0 


-1.66 


-6.00 


-6399.332 


-.0243 


135.457 


-34632 


1030.0 


-3.88 


-1.26 


-23040.770 


-.0063 


107.509 


.060164 


1339.0 


-2.22 


-0.43 


-64212.942 


-.0022 


171.079 


.709676 


1081.5 


-2.77 


-0.82 


-34474.195 


-.0041 


111.440 


.01238 


1442.0 


-2.77 


-0.46 


-59822.933 


-.0024 


146.968 


.145697 


515.0 


-1.11 


-1.44 


-20410.493 


-.0071 


54.025 


.126881 


1442.0 


-6.11 


-1.01 


-28242.996 


-.0050 


149.440 


.06137 


1133.0 


-2.77 


-0.75 


-37639.751 


-.0038 


116.493 


.042618 


618.0 


-2.22 


-1.98 


-15405.752 


-.0095 


65.760 


.095000 


1133.0 


-l.ii 


-0.29 


-92272.089 


-.0015 


114.868 


.111753 


2060.0 


-2.22 


-0.18 


-150915.860. 


- . 0009 

1 


208.067' 


.158484 


875.5 


-4.44 


-1.98 


-15426.860 


-.0095 


93.154 


.089827 


772.5 


- .55 


-0.32 


-85614.909 


-.0017 


78.378 


.098703 


618.0 


-1.11 


-0.99 


-28764.428 


- . 0049 


64.009 


.035126 


721.0 


-3.88 


-2.55 


-12423.038 


-.0119 


77.801 


.208904 


875.5 


-1.66 


-0.74 


-37875.237 


-.0038 


55.825 


.174716 


1133.0 


-2.79 


-0.74 


-37980.775 


-.0038 


116.468 


.098703 


1030.0 


-2.22 


-0.71 


-39240.355 


-.0036 


105.800 


.098703 


721.0 


-3.33 


-2.18 


-14163.219 


-.0103 


98.529 


.738132 


1133.0 


-0.56 


-0.15 


-184229.295 


-.0008 


104.304 


.062660 


1133.0 


-3.89 


-1.02 


-29043.613 


-.0051 


117.4444 


.098703 


1287.5 


-2.78 


-0.56 


-49075.979 


-.0029 


147.914 


.315771 


1081,5 


-1.11 


-0.32 


-85030.331 


-.0017 


109.737 


.018613 


1236.0 


-1.11 


-0.25 


-110239.246 


=-0013 


125.113 


.057789 


1545.0 


-6.11 


-0.86 


-32759.215 


- . 0(344 


148.858 


.087884 


772.5 


-1.11 


-0.63 


-44052.566 


- „ 0032 


79.153 


.098703 


618.0 


-1. 67 


-1.48 


-19944.649 


-.0073 


37.485 


.023398 


721.0 


-3.89 


-2.54 


-12462.568 


-.0118 


65.775 


.103646 
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